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Abstract 

The problem of <i-dimensional Schrodinger equations with a position-dependent 
mass is analyzed in the framework of first-order intertwining operators. With the 
pair (H, Hi ) of intertwined Hamiltonians one can associate another pair of second- 
order partial differential operators (R, R\), related to the same intertwining operator 
and such that H (resp. H\) commutes with R (resp. R\). This property is inter- 
preted in superalgebraic terms in the context of supersymmetric quantum mechanics 
(SUSYQM). In the two-dimensional solution to the resulting system of partial 

differential equations is obtained and used to build a physically-relevant model de- 
picting a particle moving in a semi-infinite layer. Such a model is solved by employing 
either the commutativity of H with some second-order partial differential operator 
L and the resulting separability of the Schrodinger equation or that of H and R to- 
gether with SUSYQM and shape-invariance techniques. The relation between both 
approaches is also studied. 
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1 INTRODUCTION 



The concept of position-dependent mass (PDM) is known to play an important role in the 
energy- density functional approach to the quantum many-body problem in the context of 
nonlocal terms of the accompanying potential. This formalism has been extensively used 
in nuclei [T], quantum liquids [2|, 3 He clusters j3J, and metal clusters [I]. 

Another area wherein the PDM approximation provides a very useful tool is the study 
of electronic properties of many condensed-matter systems, such as semiconductors jS] and 
quantum dots |§J. In particular, recent progress in crystal-growth techniques (molecular- 
beam-epitaxy technique, for instance) for producing nonuniform semiconductor specimens, 
wherein the carrier effective mass depends on position, has considerably enhanced the in- 
terest in the theoretical description of semiconductor heterostructures. 

Furthermore, PDM presence in quantum mechanical problems may reflect some other 
unconventional effects, such as a deformation of the canonical commutation relations or a 
curvature of the underlying space It has also recently been signalled in the rapidly- 
growing field of PT-symmetric JBJ (or, more generally, pseudo-Hermitian [9 ) quantum me- 
chanics as occurring in the Hermitian Hamiltonian equivalent to the PT-symmetric cubic 
anharmonic oscillator at lowest order of perturbation theory [TTij . 

All these developments have stimulated the search for exact solutions of PDM quan- 
tum mechanical problems both in the nonrelativistic |7J [TTJ [T2J EH EH H3 CHI H3 EH EH 
EDI EH 1221 ESI Elj and relativistic [23 EEj contexts since they may provide a conceptual 
understanding of some physical phenomena, as well as a testing ground for some approxi- 
mation schemes. In the nonrelativistic case, which we are going to consider here, the gen- 
eration of PDM and potential pairs leading to exactly solvable, quasi-exactly solvable or 
conditionally-exactly solvable Schrodinger equations has been achieved by extending some 
methods known in the constant-mass case, such as point canonical transformations [2Zj, 
Lie algebraic methods j2H], as well as supersymmetric quantum mechanical (SUSYQM) 
and shape-invariance techniques [2H|. 

Most of these works have been devoted to one- dimensional systems. The aim of the 
present paper is to tackle the more difficult problem of generating exact solutions for d- 
dimensional PDM Schrodinger equations. For such a purpose, we will extend an approach 
that has proved very useful in one dimension, namely that of intertwining operators [T§] . 
related to SUSYQM methods (see also [3 H21 HS1 HH HS1 CHI CHI 1201 12H ) As a first step, 
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we plan to consider here the case of first-order intertwining operators and their SUSYQM 
interpretation. 

In Section 2, the general theory of <i-dimensional PDM Hamiltonians admitting a first- 
order intertwining operator is developed and a solution to the resulting system of partial 
differential equations is obtained in the two-dimensional case. From this solution, a new 
exactly-solvable PDM model in a semi-infinite layer is built in Section 3. In Subsections 3.1 
and 3.2, the separability of the corresponding Schrodinger equation and the existence of an 
intertwining operator are shown to provide two alternative methods for getting the bound- 
state spectrum and wavef unctions. The relation between both approaches is then reviewed 
in Subsection 3.3. Section 4 contains some final remarks. 



2 GENERAL THEORY AND SUPERSYMMETRIC 
INTERPRETATION 



On using the von Roos general two-parameter form of the PDM kinetic energy operator [371] . 
choosing units wherein h = 2mo = 1 and summing over dummy indices, the (^-dimensional 
PDM Schrodinger equation may be written as 



Here x = (07, x 2 , . . . 9j = d/dx{, % = 1, 2, . . . , d, M(x) is the dimensionless form of 
the mass function m(x) = m M(x), V(x) is the potential and a, /3, 7 are the von Roos 
ambiguity parameters, constrained by the condition a + (3 + 7 = —1. The expression of 
the kinetic energy operator in (j2.1j) has the advantage of having an inbuilt Hermiticity 
and of containing as special cases all the proposals made in the literature to cope with 
noncommutativity of the momentum and PDM operators. 

As is well known (see, e.g., [3112]), we can get rid of the ambiguity parameters in the 
kinetic energy operator by transferring them to the effective potential energy of the variable 
mass system. Equation (j2.1|) therefore acquires the form 




(2.1) 




(2.2) 



where 



V cK {x) = V{x)+ l -{f3 + l) 



AM 



[a(a + P + l) + P + l] 



(diM^diM) 



(2.3) 



M 2 



M 3 
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and A = didi denotes the d- dimensional Laplacian operator. 

In the following, we shall be interested in the bound-state energy spectrum E n , n = 0, 
1, 2, . . . , of H and in the corresponding wavefunctions i/) n (x), n = 0, 1, 2, . . . . To be 
physically acceptable, the latter should satisfy two conditions 20/. 

(i) As in conventional quantum mechanics with a constant mass, they should be square 
integrable over the domain D on which M(x) and V(x) are defined, i.e., 

/ dx\^ n (x)\ 2 < oo. (2.4) 
Jd 

(ii) To ensure that H be Hermitian in the Hilbert space to which the ^> n (a?)'s belong, they 
should be such that 

lib (x)\ 2 

lv / v ;| ^ on the boundary of D. (2.5) 

This extra condition may have relevant effects whenever M(x) vanishes at some boundary 
point. 

Let us consider the intertwining relationship 

r)H = Htf, (2.6) 

where Hi has the same kinetic energy term as H and an associated effective potential 
Vi >e s(x). Choosing a first-order intertwining operator 

Tj = A®(x)di + B(x), (2.7) 

we are led to the restrictions 

+ djA^ = -StjAW^ (2.8) 

1 cff Vl ' eS) \ M 2 M3 J M 

+ ^ ^ ~ l 9 ^ (2 ' 9) 

A^diVes = —B(V e g - V lteS ) - ^AB + {diB)^. (2.10) 
From the intertwining relation (J2.6|) and the Hermiticity of H and Hi, it follows that 

Hrf = rfH u (2.11) 
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showing that 

R = ri\ i2i = W (2-12) 

are such that 

[H,R]=0, [H 1 ,R 1 ]=0. (2.13) 

Hence R is an integral of motion for the quantum mechanical problem described by H 
(provided it is independent of the latter). A similar situation occurs for the pair (Hx,Rx)- 
Inserting r] from (I2~7I) and r/ f = -A®di - fyA®) + B in (l2~T2l) . we get 

R = -A^A^dl - [A^idjA^) + AW{djA®)\di - A {l \diB) - B(diA {i) ) + B 2 , 
Rx = R-A®[($ j A®)-2(d i B)]. (2.14) 

The results obtained so far are amenable to a superalgebraic interpretation. Since 

r]R = R x -q, Rrf = rfR u (2.15) 

we may indeed introduce some matrix operators 

»-(?£)■ a+ K:o)' 

(2.16) 

satisfying the defining relations of the gl(l/l) superalgebra ; 2I], i.e., 

{Q + , Or} = 0, {Q ± , Q ± } = 0, [TZ, Q±] = [ft, Q±] = [ft, ft] = 0. (2.17) 

In such a context, 7?., <2 + and Qr realize the usual SUSYQM si (1/1) superalgebra [2T?] . 
while 7i turns out to be the Casimir operator of gl(l/l). 

In the one-dimensional case, it has been shown elsewhere [HI! that the solution of 
Eqs. (JUJ) - (M can be written as 

A(x) = M~ 1 /\ V eS (x) = -(AB)' + B 2 + A, Vi )eff (x) = V eS - A(A" - 2B'), (2.18) 

where A is some integration constant. Note that from ()2.18|) onwards, we denote by a prime 
the derivative of any function depending on a single variable. On combining now (|2.18|) 
with (|2.14p . we arrive at the trivial outcome 

R = H-X, Rx = H x -X. (2.19) 
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To get significant results in (|2.12|) - (|2.17J) . we must therefore be at least in two dimen- 
sions. In such a case, the existence of the integral of motion R (resp. Ri) means that H 
(resp. Hi) is integrable. The first constraint (|2.8|) can now be written as 

d x AW = -^(AV>d x M + A®d v M), d v AW = -d x A®, d x A™ = d y A^\ (2.20) 

where we henceforth denote xi, X2, d\, 82 by x, y, d x , d y , respectively. Conditions ()2.20|) 
imply that 

AA {1) = AA {2) = 0, (2.21) 

thus leading to a small simplification in the second constraint (J2.9|) . 

Let us find the general solution to Eq. ()2.20|) in the special case where the PDM depends 
on a single variable, e.g., 

M = M(x). (2.22) 

The first constraint in ()2.2()j) then amounts to 

d x A^ M' 
AW ~ ~2M' 

so that 



(2.23) 



A^ = (2.24) 
y/M 



where f(y) is a so far undetermined function of y. Plugging (J2.24j) in the second constraint 
of fl£HJ), we find 

in terms of another undetermined function g(y). It remains to solve the third constraint in 
(|2.20j) . which can be rewritten as 

On deriving both sides of this relation with respect to x, we arrive at the conclusion that 

m (wggT _ .„_. 

where C is some constant. Let us now distinguish between the three subcases C < 0, C = 
and C > 0. 
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Whenever C = —q 2 < (where we may assume q > 0), we find from ()2.27|) 

f(y) = asinqy + b cos qy, — = = csinh gx + dcoshqx, (2.28) 



where a, b, c, d are four integration constants. Inserting these relations in ()2.26|) shows that 
g(y) must reduce to some constant g. We conclude that in the first subcase, the general 
solution of (|2.2(Jj) . corresponding to the choice (|2.22|) . is given by 

^ ^ (csinh qx + dcoshqx) 2 ' 
A^\x,y) = (asinqy + b cos qy) (csinh qx + d cosh qx), 

A^ 2 \x,y) = (— a cos qy + b sin qy)(c cosh qx + dsinhqx) + g. (2.29) 

By proceeding in a similar way for C = and C = q 2 (with q > 0), we respectively 
obtain 

M(x) -- ' 



(cx + d) 2 ' 
A (1) (x,y) = (ay + b)(cx + d), 

A {2) (x, y) = --ac(x 2 - y 2 ) - adx + bey + g, (2.30) 



and 



M(x) 



(c sin qx + d cos qx) 2 ' 
A^\x,y) = (a smb. qy + b cosh qy)(c sin qx + dcosqx), 

A^ 2 \x,y) = (a cosh qy + bsinhqy)(ccosqx — dsinqx) + g, (2-31) 

in terms of five arbitrary constants a, b, c, d, g. 

Considering next Eq. (j2.9J) in the case where M is given by (|2.22|) . we get 

1C 2 

A^Ote-v^) = -(±Xa.A<»-^ayB, (2.32) 

where use has been made of Eq. (12.27)) . 

To go forward in the solution of such relations, it is appropriate to make some specific 
choice for C and for the five integration constants a, b, c, d, g. Here we consider the first 
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subcase referred to hereabove and therefore assume C = —q 2 < 0. Some comments on 
the two remaining subcases will be postponed until Section 4. On taking a — d — 1 and 
b = c = g = 0m 1)2.29)1 , we arrive at 

M(x) = sech 2 gx, A^\x, y) = cosh qx sin qy, A^ 2 \x,y) = — smb. qx cos qy. (2.33) 

Note that with this choice, the mass is well defined on the whole real line (which would not 
happen if c did not vanish). Furthermore, q may be treated as a deformation parameter so 
that when q — > 0, we get back the constant-mass case M(x) — > 1. 
The two conditions in (j2.32j) can now be rewritten as 

Kff — ^i,eff — 2g 2 cosh 2 qx — 2 cosh qx esc qy d x B (2.34) 

and 

V c s — Vi ofj = — 2g 2 cosh 2 qx + 2 cosh qx coth qx sec qy d y B, (2.35) 

respectively. On equating their right-hand sides and setting B(x, y) = sinqy (3(x, y), we get 
a separable first-order partial differential equation for /3(x,y). Its general solution leads to 

B(x, y) = (q sinh qx + F csch qx) sin qy + G, (2.36) 

where F and G denote some integration constants. From (|2.34|) and ([2.36)1 . we then obtain 

V e g - Vijtf = 2qF coth 2 qx. (2.37) 

Finally it remains to solve Eq. (j2.10jl . On taking ()2.33jh ()2.36)) and ()2.37)) into account, 
it can be easily transformed into the partial differential equation 

coth qx d x V cS - cot qy d y V eS 
= —2q[q 2 cosh 2 qx + F(q + F) csch 2 qx coth 2 qx + FG csch qx coth 2 qx esc qy] .(2.38) 

We observe that such an equation is separable provided we choose a vanishing integration 
constant G. In this case, we obtain 

V e g(x, y) = —q 2 cosh 2 qx + F(q + F) csch 2 qx + Jln(coshga; sec qy) + K, (2.39) 

Vi,eft( x , y) — ~Q 2 cosh 2 qx + F(—q + F) csch 2 qx + J ln(cosh qx sec qy) + K — 2qF, (2.40) 
in terms of two additional integration constants J and K. 
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We conclude that Eqs. lj£3fil) (with G = 0), (jQH) and (g3D|) provide us with a 

particular solution to the system of partial differential equations ()2.8j) - ()2.f 0)1 in the d = 2 
case. In order that this solution be viable and useful in the context of PDM Schrodinger 
equations with a bound-state spectrum, we have to impose some additional restrictions. 
This will be the purpose of Section 3, where we shall proceed to construct and to solve an 
interesting two-dimensional PDM model. 

3 A NEW EXACTLY-SOLVABLE PDM MODEL IN 
A SEMI-INFINITE LAYER 

Let us consider a particle moving in a semi-infinite layer of width n/q, parallel to the x-axis 
and with impenetrable barriers at the boundaries. This means that the variables x, y vary 
in the domain 

7T 7T 

D : < x < oo, < y < — , (3.1) 

2q 2q 

and that the wavefunctions have to satisfy the conditions 

V(0,y) = 0, VU±|J=0. (3.2) 

The motion of the particle is assumed to be described by the Hamiltonian 

H^ = —d x cosh 2 qx d x — cosh 2 qx d 2 — q 2 cosh 2 qx + q 2 k(k — 1) csch 2 qx + q 2 Vo, (3.3) 

corresponding to the mass and the effective potential given in (j2.33J) and ()2.39|) . where in 
the latter equation we have set J = to get an acceptable potential, as well as F = —qk 
(with k > 0) and K = q 2 Vo (with vq arbitrary). 

From the results of Section 2, it follows that and 

Hi = —d x cosh 2 qx d x — cosh 2 qx d 2 — q 2 cosh 2 qx + q 2 k(k + 1) csch 2 qx + q 2 (v + 2k) (3.4) 

admit the intertwining operator 

rf k > = sin gy(cosh qxd x + q sinh qx — qk csch qx) — sinh qx cos qy d y (3.5) 

and that their respective constants of motion read 

= — cosh 2 qx sin 2 qyd 2 + 2 sinh qx cosh qx sin qy cos qy d 2 y — sinh 2 qx cos 2 qy d 2 
+ q sinh qx cosh qx(l — 4 sin 2 qy)d x + q sin qy cos qy(l + 4 sinh 2 qx)d y 
+ g 2 (sinh 2 qx — sin 2 qy — 3 sinh 2 gx sin 2 qy) — q 2 k(l + csch 2 sin 2 gy) 
+ g 2 /c 2 csch 2 qx sin 2 (3.6) 
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and 

R (k) = R (k) + 2( fk{l + csch 2 qx sin 2 qy). (3.7) 

The bound-state energy spectrum and wavefunctions of can be determined in two 
alternative ways: either by directly solving the corresponding Schrodinger equation, which 
is separable in the variables x, y, or by taking advantage of the intertwining operator rj^ k \ 
We plan to successively review both approaches. Then, in a third step, we will determine 
the relation between the two resulting sets of wavefunctions. 

3.1 Separation of Variables in the PDM Schrodinger Equation 

Let us consider the Schrodinger equation (|2.2j) for with corresponding energy eigenval- 
ues and wavefunctions denoted by and tjj^ k \x,y), respectively. From ()3.3|) . it is clear 
that the operator 

L = -d 2 (3.8) 

commutes with H^ k \ so that both operators are simultaneously diagonalizable or, in other 
words, Eq. ()2.2|) is separable in the variables x, y. 

The dependence of the wavefunctions on y is determined by the solution of the eigenvalue 
problem for L on the interval —7t/(2q) < y < ir/(2q) with boundary conditions coming from 
the second relation in (J3.2j) . One easily finds 

L X i(y) = (l + 1) V*i(?)> 1 = 0,1,2,..., (3.9) 
where the (normalized) eigenfunctions are given by 

f J%coa[(! + l)qy] for I = 0, 2, 4, . . . 

X i(y) = < V • ( 3 - 10 ) 

[ y/%sm[(l + l)qy] for I = 1, 3, 5, . . . 

The simultaneous eigenfunctions of and L can therefore be written as 

^,y) = ^{x)xi{y), (s.ii) 

where 4>^}(x) satisfies the ordinary differential equation 

d 2 d 
cosh 2 qx—^r — 2q sinh qx cosh qx— — h q 2 l(l + 2) cosh 2 qx + q 2 k(k — 1) csch 2 qx 
dx z dx 

{e^J - q 2 v ) ^}(x) (3.12) 
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and n distinguishes between the independent solutions characterized by the same eigenvalue 
of L. 

By the changes of variable and of function 



qx 



tan 



1 / 7T 

- [z + - 

2 V 2 



71 

< z < -, 
2' 



(3.13) 



4> { $[x{z)]=VSS^&> } *(z), (3.14) 

Eq. (J3.12)) can be mapped onto the (constant-mass) Schrodinger equation for the two- 
parameter trigonometric Poschl- Teller potential (also termed Poschl- Teller I potential) [32] 

d 2 



where 



dz 2 



+ k{k - 1) esc 2 z + A(A - 1) sec 2 z &' x \z) = {z) 



(3.15) 



n ~ q 2 



E 



(k) 



q 2 v + q 2 [k- - 



(3.16) 



From the known solutions of the latter [2H1 E3J in terms of Jacobi polynomials, we obtain 
for the former 

E ni = q 2 K 2n + l + 2)(2n + l + 2k + l)+ v ] 



(3.17) 



and 



Af${t&nhqx) k {sechqx) l+2 P r [ k - 2tanh 2 gx), 

2g(2n + j + fc + |)n! r(n + Z + k + |) x 1/2 
(n + / + l)!r(n + fc + i) 



(3.18) 



The functions f|3. 1 lj) with (j)^\(x) and X/(y) given in Eqs. ()3.18|) and ()3.10|) . respectively, 
are square integrable over D with the normalization factors as indicated. Since the PDM 
M{x) = sech 2 qx vanishes for x — *> oo, such functions will describe bound states only if 
condition ()2.5|) is satisfied in this limit. This is indeed the case because coshgx| 



x ~ 



exp[— (2/ + 3)qx] for x — > oo. 

We conclude that the bound-state spectrum of is made of an infinite number of 
energy levels, given by (|3.17|) . Since E^} only depends on the combination 



N = 2n + /, 



(3.19) 
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we can rewrite it as 

= E$ = q 2 [{N + 2){N + 2k + 1) + v ]. (3.20) 

The level characterized by iV comprises the states specified by (n, I) = (0, N), (1, N—2), . . . , 
(y, 0) or (^p, 1) according to whether N is even or odd. Its degeneracy is therefore given 
by 



deg(iV) 



N 
~2~ 



+ i, (3.21; 



where [iV/2] stands for the integer part of N/2. 

Before going to the alternative approach based upon the intertwining operators r/^ and 
^( fc )t m Subsection 3.2, it is worth making a pause for observing that the one-dimensional 
Hamiltonian H\ , as defined in (j3.12|) . itself satisfies some intertwining relationships 

A (k) H (k) = r H (k+i) + 2q 2j\ A (k)^ A (k) H ( k ) = r H (k+i) + 2q2} \ A {k) (3 22) 

with the first-order differential operators 

A^ = cosh qx— — h q(l + 2) sinh qx — qk csch qx, (3.23) 
ax 

A^ = cosh qx ql sinh qx — qk csch qx. (3.24) 

ax 

(k) 

Such a property derives from the existence for H t of a four-way factorization 

rrM _ A {k)\ Ak) (k) _ ,(k-l) .(fc-l)f , (fc) 

_ 7(*)t J(*) , _ I(*-l)t , 

c[ fc) ee g 2 (Z + 2)(/ + 2A; + l), c[ fe) ee g 2 Z(Z - 2A; + 1). (3.25) 

The latter is similar to that occurring for the three-dimensional radial harmonic oscilla- 
tor [34J and is directly connected with the existence of an so(4) potential algebra for the 
two-parameter trigonometric Poschl- Teller potential [3*3*] . 

In a standard (unbroken) SUSYQM approach to H^, one would only consider the pair 
of operators A\ and A\ ; . The former indeed annihilates the ground-state wavefunction 
0Q fc / of and changes cp^\ into <j>n-i i+\ for n = 1, 2, . . . , while the latter generates 
transitions between anc ^ ^i+i j for n = 0, 1, 2, . . . . In contrast, the second pair of 

operators Ai and A-i does not change the n value as one goes this time from (jy*\ to 
<fini~-l for = 0, 1, 2, . . . , or vice versa. 
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Interpreting these results in terms of the full two-dimensional Hamiltonian is dif- 
ficult because one is confronted with the same type of problem as that occurring in the 
SUSYQM approach to the Coulomb problem [2H]. In such a case, the intertwining opera- 
tors for the radial Hamiltonian change the angular momentum quantum number / by one 
unit without affecting the angular part Yi m (9,(p) of the wavefunction. As a consequence, 
there has been much controversy on how to apply the SUSYQM results to atomic spec- 
tra [35J. Here too we may observe that both intertwining operators A[ k \ JO} and their 
Hermitian conjugates change by one unit the quantum number /, now connected with the 
eigenvalue of the operator L defined in ()3.8)L although the corresponding eigenfunction 
Xi(y) remains the same. 

As we plan to show in Subsection 3.3, the intertwining operators 77 w and offer the 
advantage of being free from such a problem since they act on Xi(y) as wen as 011 ( t ) ni( x )- 

3.2 Intertwining Operator Approach 

As established in Section 2, the operators H^> and defined in (13 .3|) and (|3.6j) . respec- 
tively, are simultaneously diagonalizable. Let us denote their simultaneous eigenfunctions 
by ^f^' No (x : y) and the corresponding eigenvalues by and r„ , with iV = iV + v. Hence 

H (k) ^ ] , No (x, y) = EP^ No (x, y), R^P, No (x, y) = rf)*<\(x, y), (3.26) 

where r„ > since is a positive-definite operator. We plan to show herebelow that 
in ()3.26|) . N and v run over all nonnegative integers, while N Q is restricted to nonnegative 
even integers. 

Our first step consists in constructing the functions jv ( x > v) corresponding to E^' Q 
and = 0. Since they belong to the subspace spanned by the zero modes of r/^, let us 
inquire into such zero modes and therefore consider the equation 

i ? «w?)(ar,y) = 0, (3.27) 

where s will serve to distinguish between independent solutions. From ([3.5)1 . it follows 
that this first-order partial differential equation is separable. Up to some multiplicative 
constant, its general solution is given by 

ujf\x,y) = (taioh qx) h (sech qx) s+1 (cos qy) s , (3.28) 



13 



where s is related to the separation constant. To get functions satisfying the boundary 
conditions (|3.2|) . we have to restrict s to positive values. It is then easy to check that 
u>i k \x,y) is normalizable on the domain D, defined in (13. ip . and satisfies condition (|2.5|) 
for x — ► oo. 

On applying on the zero modes (|3.28|) . we obtain after a straightforward calculation 

H^ k \x, y) = q 2 [{s + l)(2k + s) + i*]w«(a;, y) - q 2 s(s - l)coi%(x, y). (3.29) 

Hence, amongst the zero modes u)s (x,y) with s > 0, there is only one that is an eigen- 
function of H^ k \ namely u[ k \x,y): 

H^J 1 k \x,y)=q 2 [2(2k+l) + v ]J 1 k \x,y). (3.30) 

As it can be checked by comparison with (|3.20|) and with (j3.10|) . f!3.11|) . ()3.18|h the eigenvalue 
and eigenfunction of in ()3.30|) coincide with the ground-state energy Eq and ground- 
state wavefunction ip^l(x,y), respectively. We have therefore found the (unique) solution 
of (HHj) with r {k) = and N = N = 0, 

V$(x,y) =4%,V) <xui k \x,y). (3.31) 

We can now get the remaining solutions of (m with r {k) = and iV = N ^ by 
looking for those linear combinations of ui (x, y), s = 1, 3, . . . , N + 1, 

Nq+1 

*%U(*>V) = E ot 1 " (3-32) 



2' 

8=1 



that are eigenf unctions of H^ k \ Note that in ()3.32|) . N is necessarily restricted to even 
integers. On using (|3.29p . we obtain a recursion relation for the coefficients a« together 
with the eigenvalue of H^ k \ 

E { nI = Q 2 K N o + 2)(iVo + 2k + 1) + v ). (3.33) 

The solution to the recursion relation reads 
flW _ , l)( Ar 0+ l-8)/ 2 (iVo + l)!r(^ + A: + l) 

(3.34) 

The remaining multiplicative factor a N +1 can in principle be determined through the nor- 
malization condition of fy^ oNo (x,y) on D. However, since the zero modes ui (x,y), s — 1, 
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3,5,..., form a nonorthogonal set, this condition is not easy to work out. We shall therefore 
leave the construction of fully normalized functions No (x,y) to Subsection 3.3, where 
they will be expressed in an orthogonal basis. 

In a second step, we can now obtain from the functions tyj^ N (x, y) the solutions of 
(j3.2(ij) associated to v = 1, 2, . . . , and correspondingly N = AT + 1, Nq + 2, . . . , by a 
straightforward application of the intertwining relations ()2.6|) and f)2.15j) for the operators 
v (h) } r/ (fc)t ; jj(k) t H V») t R (k) t and R (k) on comparing (jS3) and flUj) with (Q and (jHSJ), 
respectively, we indeed observe that 

H (k) = H (k+i) + e ( fc)j £ (*) = 3^ (3>35) 

In other words, iif^ and i?^ are both shape invariant [2H|- in view of this property, we 
directly obtain 

v) = A4V (fe) V fc+1)t • • • ^^ffi&fo y). ^ = i,2,..., (s.sr) 

where A/jy ]y ^ s some normalization coefficient and 

Ej, fc) = ^ o +i/) + ^e (fc+l) = g 2 [(iV + 2)(iV + 2A; + l)+ Vo ], (3.38) 

i=0 

rf) = = g 2 z/(z/ + 2Jfe), (3.39) 

i=0 

for N — Nq + u. Note that the eigenvalues ()3.38|) of coincide with those previously 
found in ()3.20|) as it should be. However, the supersymmetric approach adopted in this 
Subsection has the advantage of explaining the dependence of these eigenvalues on a single 
quantum number N, which otherwise would appear to be accidental. 

As a final point, let us observe that the normalization coefficient A7^]y in (|3.37|) can be 
easily derived from the property 

/^(*) ^(*) \ _ [ ^ /^(fc+i) |„(*)„(fc)t I \ f3 4n) 

\ JV N-1,N / 

and Eqs. (l2~T2l) . (I3~2T)J) . (JSSi, B- The result reads 
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3.3 Relation between the Two Approaches 

In Subsections 3.1 and 3.2, we followed two different methods for determining the bound- 
state spectrum and wavefunctions of H^ k \ The former was based on the separability 
and exact solvability of the corresponding Schrodinger equation, while the latter used an 
intertwining-operator approach and its supersymmetric interpretation. Such methods re- 
sulted in two distinct bases for the associated Hilbert space, ip^j(x,y), n, I = 0, 1, 2, . . . , 
and $?P jNq (x, y), N = 0, 2, 4, . . . , N = N , N + 1, N + 2, . . . , corresponding to definite 
eigenvalues of the respective symmetry operator, L or R^ k \ The aim of this Subsection is 
to determine the transformation matrix between these two bases. 

For such a purpose, let us first consider the action of the intertwining operators 77 w 
and on the first basis wavefunctions ip^j(x,y). As it is clear from the analysis of 
Subsection 3.2 (see Eq. ()3.37|) ). such operators have the property of decreasing or increasing 
iV = 2n+l by one unit, while changing k into k+1 or vice versa. Hence, when acting on some 
generic function ij)^\(x,y), r]^ k \ for instance, will generate transitions to ( x ' v) anc ^ 

i^ni^i( x ^y)- This is confirmed by an explicit calculation using Eqs. (|3.5|) . (|3.1U|) . (|3.11|) . 
f)3. 18j) . and some well-known properties of Jacobi polynomials [36J. The obtained result 
reads 



We are now able to expand the members of the second basis that are zero modes of rf k \ 
i.e., W$ nX x i v)i No — 0, 2, 4, ... , into linear combinations of the first basis wavefunctions 




(3.42) 



Similarly, for the Hermitian conjugate operator we get 




(3.43) 



T*j(x,y) with 2n + I = iV c 



N /2 



.N -2n 



€,N -2n(^Vl AT =0,2,4,..., 



(3.44) 



n=0 
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where X^ No _ 2n are some coefficients to be determined. On applying 77^ on both sides of 
this equation and using the facts that the left-hand side vanishes while Eq. f!3.42|) allows 
us to calculate the right-hand one, we indeed obtain a recursion relation for the coefficients 
Xn^N -2m whose solution is given by 



X W I (f)'(f + l)!r(iVo-n + A: + |)r(n + A: + |) \ 1/2 (fc) 



As shown in the Appendix (see Eq. (|A.3J) 1. the remaining multiplicative factor X^\ 20 can 

(k) / \ (k) 

be computed from the normalization condition of Nq (x, y). The final result for X^ No _ 2n 
can then be written as 

1/2 



w , (N + l)\T(k + l)T(n + k + ^T(N -n + k + l) 

A »- v "- 2 " I 2^n!(JV -n+l)!r(^ + ife + l)r(ifc + i)r(^t3+ife) 1 ' °' " 

For the first few Nq values, Eq. 1)3.44)1 reads 

<l = 4% (3-47) 



= ^==T \A + § C + W3 (* + \ ) 4 k) o ) , (3-48) 



i\ / . 3 \ (fc) 



+ \l^[k + -j (k + -)^>). (3.49) 

Note that Eq. ()3.47j) coincides with (j3.31|) as it shoud be. 

To be able to extend expansion ()3.44|) to all members ^^ No (x, y) of the second basis, 
defined in ()3.37j) . we have first to iterate Eq. (|3.43|) to determine the action of v operators 
of type T)( k " on the first basis wavefunctions: 



V 



n (*)t n (W)t . . . „(fc+"-i)t A k +») _ V Y [k) tb {k) f3 501 

'/ '/ '/ V n ,l — / j 1 n,l;n+u,l+v-2uVn+u,l+v-2u- {O.OU) 

fi=0 

Here the coefficients are given by 

{n + fi)\{n + l + u-fi + iy.r(n + u + k + \)r(n + l + u + k + l)\ 1/2 f _ 
X ' n\ {n + I + 1)! T (n + /1 + k + |) f (n + I + v - Ji + k + §) ' ' 
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where x^il^v) is a combinatorial factor giving the number of ways one can go from I to 
I + v — 2fi in v steps so that in each step I increases or decreases by one unit without 
reaching negative values. This means in particular that 

~Z + v~ 



x^(l, v) = if jj, > 



(3.52) 



Hence the summation on the right-hand side of (|3.5U|) actually runs from to min (y, [^p] ) • 
For v = 1, Eq. (|3.5U|) reduces to Eq. (|3.43|) since in such a case 

x (l,l) = l, x 1 (l,l) = l-6 lfi . (3.53) 

For higher v values, it can be easily proved by induction over v on taking the relations 

x (l,i>) = x (l,u-l), 

x^l, v) = (l - S^ ( i + „ +1) / 2 ) x^(l,u - 1) +x fl (l,v- 1), n = l,2,...,v-l, 
x v {l,v) = (l-VOx^^iz-l) (3.54) 

into account. Some examples of explicit solutions to this set of recursion relations are 
provided in the Appendix. 

Finally, it remains to combine Eqs. ()3.37|) . ()3.44|) . and ()3.50|) to obtain the searched for 
expansion 

[N/2] 

^N,N = % >N ;n,N ~-2ni 'n,N -2W (3.55) 
n=0 

where 

n' 

' "max 

^N ;n,N-2n — JV N,N Z-e n',N -2n' 1 n',N -2n';n,N-2n- {O.OV) 

Here n' min = max(0, n+No—N), n' ma:K = min nj and the three factors on the right-hand 
side are given by Eqs. ()3.41|) . (j3.46|) and ()3.51|) . respectively. 

For those wavefunctions of the second basis ty$ (x,y) that can be reached from the 
ground state ()3.47|) by applying N operators of type , Eq. ()3.56|) can be written explic- 
itly as 

Z^ N . 2n = (-l^ N{N+3)+n (N-2n + l) 

NlT^ + k)T(-±l + k)T(N + k + l) V /2 
" 1 2 N n\ (N -n + 1)! T(N + k)T (n + k + ±) T (N - n + k + |) 1 
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on using Eq. (jA.8|) for the combinatorial factor x^(0, v). So, for instance, 

^1=4% (3-58) 




r(*) - 1 I . U fix M ./.(*) , ./l. . 5 „/.W 




= ^=y ( - + 5 <3 + aA + J ^5 ) , (3.60) 



+ J 2 ( k+ l) ( k+ l)^) ■ (3 - 61) 

As it can be checked, the functions (13.48)1 and (|3.59J1 . as well as (|3.49jl and (|3.61|) . are 
orthogonal. 

4 FINAL REMARKS 

In this paper, we have developed a general framework for analyzing (^-dimensional PDM 
Hamiltonian pairs (H, Hi) admitting a first-order intertwining operator rj. We have estab- 
lished that for d > 2, there always exists another pair of intertwined second-order partial 
differential operators (R, Ri) associated with the same operator rj and such that R (resp. 
Ri) commutes with H (resp. Hi). We have shown that in a SUSYQM context based on 
an sl(l/l) superalgebra, R and R\ can be interpreted as SUSY partners, while H and Hi 
are related to the Casimir operator of a larger gl(l/l) superalgebra. Furthermore, we have 
derived a system of partial differential equations to be satisfied by the mass, the functions 
appearing in the definition of 77 and the effective potentials contained in H and H\. 

Considering this system in more detail in the two-dimensional case for a PDM depending 
on a single variable, we have proved that the mass may be of a hyperbolic, trigonomet- 
ric or rational nature. We have then chosen the class of hyperbolic PDM's and, more 
specifically, a sech 2 -mass, which has proved very useful in the context of exactly-solvable 
one-dimensional problems |T9"1 121j. Under this assumption, we have obtained the general 
solution to the system of partial differential equations when some integration constant (G 
in (|2.36jl ) vanishes. 
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From such a solution, we have built a physically-relevant model depicting the motion of a 
particle in a semi-infinite layer. The corresponding Hamiltonian H has two noncommuting 
integrals of motion, L and R. Diagonalizing H and L simultaneously amounts to separating 
the Schrodinger equation into two exactly-solvable differential equations. On proceeding 
along these lines, we have obtained a bound-state spectrum made of an infinite number 
of levels exhibiting some finite degeneracies. Choosing next to simultanously diagonalize 
the shape-invariant operators H and R, we have taken advantage of poweful SUSYQM 
techniques to obtain their respective spectrum and eigenfunctions in a straightforward way. 
The latter approach has had two interesting outcomes. First, it has shed some light on 
the origin of the 'accidental' degeneracies observed in the Hamiltonian spectrum derived in 
the former approach. Second, it has provided us with a partnership between the full two- 
dimensional Hamiltonians H and Hi, free from the interpretation problems that plague the 
more traditional description based on one- dimensional intertwining operators [35] . 

Had we chosen the class (|2.3H1 of trigonometric PDM's instead of that of hyperbolic 
ones, a solution to the system of partial differential equations similar to that obtained in 
Section 2 would have emerged. The counterpart of the model considered in Section 3, now 
depicting a particle moving in a rectangular box, however turns out to be devoid of bound 
states. Constructing some physically-relevant model in such a case therefore remains under 
study. 

Finally, in view of an equivalence between Schrodinger equations involving a PDM and 
those in a curved space [7], selecting the class ()2.30|) of rational PDM's and, more specifically, 
a l/x 2 -mass would have led to some results directly connected with those recently obtained 
for intertwined Hamiltonians in Poincare half plane (AdS-z) [3*7j . 

From a physical viewpoint, our model might find applications in the study of quantum 
wires with an abrupt termination in an environment that can be modelled by a dependence 
of the carrier effective mass on the position. The presence of bound states may be compared 
to that observed in a quantum channel whenever the uniformity is broken, for instance by 
a quantum dot or a bend (see, e.g., jSH] and references quoted therein). 

Some interesting open problems for future work are the extensions of the present study to 
higher-dimensional models and to second-order intertwining operators (for a recent review 
on the latter see, e.g., [39 J. 
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APPENDIX 



In this Appendix, we prove some results used in the calculation of the transformation matrix 
between the two bases {ip^j} and {^^ No } in Subsection 3.3. 

Let us first consider the computation of the coefficient X^\ 2Q in expansion ()3.44j) . On 

using Eq. ()3.45|) . the normalization condition Yln=o (^njv -2n) = 1 of ^f^l,N ( x i v) can 
be written in terms of 

N ° ~ ^ n! (JVo-n + 1)! ~ ^ = n! (jV -ra + 1)! ' 

(A.l) 

The identity (1 + x) _ ( fc+ ^)(l + = (1 + x) _(2fe+1) , combined with the binomial 

theorem, leads to the closed form 



w r(iv + 2fc + 2)r 2 (fc + 1) 

*° " 2(iV + l)!r(2A; + l) ' 



From this result and some elementary properties of the gamma function , we obtain 

1/2 



(N Q + l)\T(^ + k)T(k + l) 
2 n (JSi) ! (JSi + 1)! r (Jfe + 1) r (f + A; + 1) 



-^Vq/2,0 ~~ ( (V„ / Nn \ I /TVn , i^if/i. 1\ f /JV» , L , 1\ I ' (A-3) 



Let us next provide some solutions to the set ()3.54j) of recursion relations for the com- 
binatorial factor Xfj,(l, v) with the starting values given in ()3.53j) . The simplest case cor- 
responds to I > v, because in none of the v steps from I to I + v — 2fi can a negative I 
value be reached. From (|3.52j) . it follows that x^(l, v) is nonvanishing for // = 0, 1, . . . , v. 
Furthermore, Eq. (|3.54j) reduces to 

x (l,u) = x (l,v-l), 

x^l, v) = x^-i^,!/- 1) +x fi (l,u- 1), /jl = 1,2,..., v- 1, 

x v {l,v) = x^&v-l), (A.4) 

so that Xfj,(l, v) is simply given by a binomial coefficient 

Xn{l,v) = I V \ , /i = 0,l,...,u, i£l>u. (A.5) 
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In contrast, the value / = maximizes the number of possible encounters with negative 
/ values during the v steps. The nonvanishing values of x^(l, v) now correspond to fj, — 0, 
1, . . . , [|] and satisfy the relations 



v 



x (0,u) = x o (0,u-l) 
x M (0,i/) = x M _i(0,i/- 
x 5 (0,i/) = a; w (0,i/-l), (A.6) 



x^(0,v) = x M _i(0,i/- 1) +x li (0,v- 1), = 1,2, 1, 



or 



rr (0, i/) = x (0,i/-l), 

^(0,1/) = ^_i(0,i/-l)+ ^(0,^-1), // = 1,2,...,^^, (A.7) 
for even and odd values of u, respectively. In both cases, the solution reads 

^ I [|] +!,[*] +2,. 
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